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ABSTRACT 

Active galactic nuclei (AGN) display an extreme range in the narrow emission line 
equivalent widths. Specifically, in the PG quasar sample the equivalent width of the 
narrow [Oiii] A5007 line has a range of > 300 (< O.SA to 157A), while the broad H/3 
line, for example, has a range of ten only (23A to 230A) . The strength of [O ill] A5007 is 
modulated by the covering factor (CF) of the narrow line region (NLR) gas, its density 
(rie), and ionization parameter ([/). To explore which of these factors produces the 
observed large range in [O ill] A5007 strength, we measure the strength of the matching 
narrow H/3 and [Oiii] A4363 lines, detected in 40 out of the 87 z < 0.5 PG quasars 
in the Boroson & Green sample. The photoionization code CLOUDY is then used to 
infer CF, rie, and U in each object, assuming a single uniform emitting zone. We find 
that the range of CF (~ 0.02 — 0.2) contributes about twice as much as the range in 
both TT-e and U towards modulating the strength of the [O ill] A5007 line. The CF is 
inversely correlated with luminosity, but it is not correlated with L/ L-^dd as previously 
speculated. The single zone [Oiii] A5007 emitting region is rather compact, having 
-Rnlr = 40L44*^ pc. These emission lines can also be fit with an extreme two zone 
model, where [Oiii] A4363 is mostly emitted by a dense {n^ — 10^ cm~'^) inner zone at 
R 



NLR 



^44^ pc, and [Oiii] A5007 by a low density (Ue = 10^ cm extended outer 



75OL44 pc. Such an extended [O ill] A5007 emission should be well 



zone at i?NL*R 

resolved by Hubble Space Telescope imaging of luminous AGN. Further constraints 
on the radial gas distribution in the NLR can be obtained from the spectral shape of 
the IR continuum emitted by the associated dust. 
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1 INTRODUCTION 

Active galactic nuclei (AGN) can show both 
broad (FWHM> 1000 km s"^) and narrow 
(FWHM< 1000 km s"^) emission lines. While both 
lines are excited by the central ionizing continuum source, 
the broad Lines arise from the immediate vicinity of the 
central massive black hole (J?<0.1 pc), and the narrow 
lines arise further out, sampling the host galaxy dynamics, 
and possibly chemical composition and gas distribution as 
well. The strength and profiles of the narrow emission lines 
were studied extensively over the past 30 years, revealing 
radial gradients in the density and ionization state of the 
gas, and some characteristic trends in the line profiles (e.g. 
Heckman et al. 1981; Osterbrock 1989; Peterson 1997; 
Vanden Berk et al. 2001). 

One remarkable property of the narrow emission lines, 
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which has not been explored yet in detail, is their extreme 
range of equivalent widths (EWs). Specifically, in the Boro- 
son & Green (1992, hereafter BG92) sample of 87 2 < 0.5 
quasars from the bright quasar survey (Schmidt & Green 
1983), the narrow [O ill] A5007 line ranges in EW from 157 A 
to undetectable (< 0.5 A), i.e. by a factor of > 300, while in 
contrast the EW of the broad H/3 line in the same sample 
ranges only over a factor of 10 (23 - 230 A). The narrow line 
EW is mainly set by two factors, the covering factor (CF) 
of the photoionized gas at the narrow line region (NLR), 
and the line emissivity (assuming no foreground absorption) . 
The CF is set by the spatial distribution of the gas at the 
NLR, and by the angular distribution of the illuminating 
ionizing radiation. The latter may be set by the contin- 
uum emission mechanism, or through obscuration further 
out (e.g. by the 'torus'). The line emissivity for coUisionaly 
excited lines, such as [O ill] A5007, is set by the gas electron 
density rie, the gas temperature, and by the gas ionization 
state. The last two properties are set by the ionizing photon 
fiux, or equivalently the ionization parameter {U = n^/us, 
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where is the density of H ionizing photons). Which of the 
above factors produces the large range in EWs of the narrow 
emission lines? 

The purpose of this paper is to answer this question for 
the [O III] A5007 hnc. Wc have chosen this hnc since it is 
generally the strongest and least blended narrow emission 
line in AGN. Our approach is to measure the strength of 
the narrow [Oiii] A4363 and H/3 lines in a large sample of 
AGN, where the [O iii] A5007 line was already measured. 
The [Oiii] A4363/[Oiii] A5007 line ratio is mainly a mea- 
sure of rie, while the H/3 line is dominated by recombina- 
tion and thus mostly sensitive to the CF of the gas. The 
combination of the three linos allows us to determine the 
CF, rie, and U of the [O III] A5007 emitting gas, and thus 
find which factor most strongly modulates the strength of 
the line. This derivation relies on the simplifying assump- 
tion that H/3, [O III] A5007, and [O in] A4363 originate in a 
single zone with homogeneous gas properties. We relax this 
assumption by also making an extreme two zone approxima- 
tion, where [O III] A5007 is mostly emitted by an outer low 
density NLR component, and [Olll] A4363 by an inner high 
density component, and again explore whether the strength 
of the [O III] A5007 line is mostly modulated by the CF or 
by the gas emissivity. 

The [O III] A5007 line is one of the main components 
contributing to the set of strongly correlated AGN emission 
properties, so called 'eigenvector 1' (hereafter EVl, BG92). 
We therefore also briefly explore whether the physical pa- 
rameters controlling the strength of [O ill] A5007 are linked 
to other EVl emission properties. The determination of U 
and He allows us to calculate the distance of the NLR from 
the central ionizing continuum source, -Rnlr, and we also 
look for relations between -Rnlr and other AGN properties. 
The paper is organized as follows, in Section 2 we provide a 
brief theoretical justification for the set of lines used in this 
analysis, and explain how the NLR parameters are derived. 
In Section 3 wc describe the sample used and the narrow line 
measurement procedure. The results are described and dis- 
cussed in Section 4, and the main conclusions are provided 
in Section 5. 



2 THEORETICAL APPROACH 

The [O III] A5007 and [O ill] A4363 hues are particularly use- 
ful for constraining T and Ue of the NLR gas which pro- 
duces the [O III] line emission. Both lines arise from the same 
element, eliminating abundance effects. Furthermore, they 
both arise from the same ionization state, eliminating ioniza- 
tion level corrections. Both lines arise through coUisional ex- 
citation, and although they correspond to similar transition 
energy, they differ in two respects. First, [O ill] A5007 arises 
from excitation to the level, located 2.51 eV above the 
ground (^Po) level, while [Olll] A4363 arises from excitation 
to the level, 5.35 eV above the ground level. Thus, their 
ratio is a useful temperature diagnostic when kT < 10 eV. In 
addition, the two transitions have different radiative decay 
rates, leading to critical densities for coUisional dexcitation 
of 7.0 X 10^ cm~^ for [O ill] A5007 and 3.3 x lO'' cnr^ for 
[Oiii] A4363. Thus, the [Oiii] A4363/[Oiii] A5007 line ratio 
is also a useful density diagnostic, when the density is not 
much above or below the two critical densities (e.g. Oster- 



brock 1989). An additional advantage of the [O iii] A5007 line 
is that it is typically the main coolant of the NLR (amount- 
ing to as much as ~ 50 per cent of the total cooling), and 
thus its absolute intensity is only a weak function of the O 
abundance (through the metalicity dependence of the gas 
temperature) . 

Figure 1 presents curves of constant 
log [Oiii] A4363/[Oiii] A5007 photon flux ratio as a 
function of the electron density and temperature. The cal- 
culations arc based on a 6-levels [O ill] model (excitations 
to levels above the 6th level, ^So, do not contribute signifi- 
cantly to the two [O III] lines at the relevant temperatures) . 
The level energies were obtained from Moore (1993, p. 267), 
and radiative transition probabilities and mean collision 
strengths from Pradhan & Peng (1995, table 1 there, for 
T = K, relevant here) ^. The figure demonstrates the 
weaker dependence of the line ratio on rie at the lowest Tie , 
and weaker dependence on T as the highest T. 

In photoionizcd gas it is more natural to use U, which 
is often a free parameter in pliotoioiiization calculations, 
rather than T which is set by U and rie through detailed 
photoionization calculations. Fig. 1 presents the dependence 
of T on Tie and U, as deduced from the photoionization code 
CLOUDY (Ferland et al. 1998). It demonstrates the weak 
dependence of T on rie, at a given U, and the overall small 
range m T (1-2.5 X 10" K) when U varies from 10""* to 
10~^. Assuming a uniform Ue gas, the measured value of 
[O III] A4363/[Olll] A5007 (corresponding to a certain curve 
in Fig .1) yields the allowed range of 

Tie ^ 

function of the 

range in U. E.g., if the measured line ratio is 0.1, then the 
implied density is Ue = W^ "^^ to 10®^ cm-^ for U = 10"^ 
to 10"". As can be seen in Fig. 1, the allowed range in rie 
increases as the line ratio gets smaller. 

Breaking the rie, U, degeneracy requires an additional 
line. The optimal line for that purpose is H/3. It is a nearly 
pure recombination line at the expected NLR conditions 
(where coUisional and optical depth effects should be small) , 
and thus its intensity provides an estimate of the fraction 
of ionizing radiation intercepted by the NLR, i.e. the NLR 
CF. The [O III] A5007/H/3 ratio provides the additional con- 
straint on rie and U, required to break the degeneracy. 
The H/3 EW, together with the EW of [Olll] A5007, and 
[Oiii] A4363, allow us to determine rie, U, and the CF of 
the [O III] emitting gas in the NLR. From the observational 
aspect the H/3 line is also optimal, as it is located close to the 
two [O III] lines, allowing all three lines to be obtained in a 
single spectrum, minimizing systematic measurement errors. 
Also, the close wavelength proximity of the lines minimizes 
dust reddening effects on the measured line ratios. 

As noted above, we use the photoionization code 
CLOUDY to calculate the strength of [Olll] A5007, 
[Oiii] A4363, and H/3, as a function of Ue (from 10^ cm~* 
to 10* cm-3), U (from 10"" to 10"^), and the CF. The 
additional input model parameters for CLOUDY are a slab 
geometry with a total column of 10^^ cm~^, ISM abundances 
with grains, and a standard Mathews & Ferland (1987; her- 
after MF) AGN continuum with a break at 1 /xm (Ferland 



1 Two collision strengths, n(2p3 5g0_ ISq) and 

Cl(2p^ ^Sj, 2p^ ^02), were not available, and we therefore 
assume are equal to 1. 
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et al. 1998). The results are not sensitive to the assumed 
column since the dust opacity strongly suppresses line emis- 
sion from columns beyond 10^^ cm~^ (Laor & Draine 1993, 
hereafter LD). The effect of the ionizing continuum shape is 
explored by calculating the model results again with a dif- 
ferent ionizing spectrum, as further described in Section 4. 

Figure 2 presents the calculated line flux ratios and EW 
as a function of rie for different values of U. The upper 
panel shows the [O ill] A5007/H/3 flux ratio. This ratio de- 
creases with increasing rie due to coUisional suppression of 
the [O III] A5007 line at rie > 10^■*^ The ratio increases with 
U at low U, as the fractional abundance of O ill increases, 
and starts decreasing at logC/ > —1.5, as O gets ionized be- 
yond O III (the dust opacity suppresses emission from deeper 
layers where O is less ionized). The middle panel shows the 
[Oiii] A4363/[Oiii] A5007 flux ratio. This ratio is mostly set 
by Tie due to the different ncrit of the two lines, eis discussed 
above. The weak dependence on U results from the weak 
dependence of T on ?7 (Fig. 1). The lower panel shows the 
H/3 EW for a CF = 1. At low U the H/3 EW is independent 
of He and U, as expected for a pure recombination line. At 
log U > —2.5 the H/3 EW starts decreasing due to the in- 
creasing suppression of the line emission by absorption in 
dust within the ionized gas (Voit 1992; LD; Netzer & Laor 
1993; Ferguson et al. 1997; Dopita et al. 2002). 

Figure 3 presents theoretical curves of the 
[O III] A5007/H/3 flux ratio as a function of the 
[Oiii] A4363/[Oiii] A5007 flux ratio for different val- 
ues of Ue and U. This plot allows one to transform the two 
measured line ratios in a given object to the corresponding 
values of Ue and U in that object. The third parameter, 
the CF, is then set by the EW of either one of the three 
lines. There is a one to one transformation from the 
[Oiii] A5007/H/3 versus [Oiii] A4363/[Oiii] A5007 coordi- 
nate to the Tie versus U coordinate, as long as log U < — 2.2. 
Above that value there are two sets of solutions, one having 
log [/ < —1.5, and the other log (7 > —1.5, for a given 
set of [Oiii] A5007/H/3 and [Oiii] A4363/[Oiii] A5007 
values (see Fig. 3). The presence of two solutions results 
from the non-monotonic dependence of [O ill] A5007/H/3 
on U, mentioned above (Fig. 2, upper panel). The two 
solutions correspond to different CF values (the larger U 
requires a larger CF to produce the observed EW, 
Fig. 2 lower panel). Obtaining a unique solution requires 
additional constraints on U or CF. As further described 
in Section 3, almost none of our objects happens to lie 
within the degenerate solution regime. However, note that 
photoionization models with logZJ > —1, not explored 
here, produce lower [O ill] A5007/H,5 ratios, which will be 
degenerate with some log (7 < — 2.2 solutions. Such high 
U models require a high CF, and they imply strong high 
ionization narrow lines, which we do not study here. 

The above analysis is within the single zone approxima- 
tion, i.e. all three lines arc assumed to originate in gas with 
a given set of rie, U and CF. In reality, there is likely to ex- 
ist a distribution of values for these three parameters in the 
NLR, in which case the relative line strength will be a func- 
tion of position in the NLR (e.g. Ferguson et al. 1997). The 
assumption of a single zone with uniform gas properties may 
not yield representative mean values in this case. To explore 
the effects of a non uniform cloud population, we repeated 
the analysis assuming two populations of clouds with very 
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Figure 1. Curves of constant log [Oiii] A4363/[Oni] A5007 pho- 
ton flux ratio as a function of log rie and log T for a 6-levels coUi- 
sionaly excited Oio. The thick lines represent T{ne) for different 
values of U, as calculated by the photoionization code CLOUDY. 
Note the weaJj dependence of T on rie at a given U. 

different [O ill] A5007 and [O ill] A4363 emissivities. The first 
population is assumed to have a density Ue = 10^ cm~^, and 
thus produces a high [Oiii] A4363/[Oiii] A5007 ratio, and 
the second population is assumed to have ne = 10^ cm~^, 
and thus produces a ratio lower by a factor of ~ 100 (Fig. 3). 
Most of the [O III] A4363 emission is expected to come from 
the first population of clouds, and most [O III] A5007 from 
the second. Such a "two zone model" has six free parame- 
ters (rie, and CF for each zone), and three need to be 
set since we only have three constraints (the EWs of H/3, 
[Oiii] A5007, and [Oiii] A4363). The additional parameter 
we fix is log U = —1 for the Ue = lO' cm~^ zone, which im- 
plies this is an inner very compact zone. This rather high U 
is motivated by the mean spectral energy distribution (SED) 
of AGN. which indicates a significant covering factor of very 
hot dusty gas with a high ionization parameter. The dust 
temperature indicates it is located outside the BLR, but on 
scales smaller than the "typical" NLR. The lack of strong 
line emission from the associated gas implies U S> 10""^, to 
allow significant suppression of the line emission by the dust 
(Section 2). The remaining free parameters of the two zone 
model are the ionization parameter of the outer zone, (/out , 
and the covering factors of the inner and outer zones, CFin, 
CFout, these are set as described below in Section 3.4. 



3 DATA ANALYSIS 
3.1 The data set 

For the analysis we use the BG92 sample mentioned above 
(Section 1). This sample extends in luminosity from Seyfert 
galaxies with vLv = 3.3 x 10*^ erg (calculated at rest 
frame 3000 A using the continuum fluxes in Neugebauer et 
al. 1987, assuming Ho = 80 km s"^ Mpc"\ VLo = 1.0), to 
luminous quasars at vL^ = 1.4x lO'*^ erg . This is a com- 
plete and well defined sample, selected based on (blue) color 
and (point like) morphology, independently of the emission 
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Figure 2. The solid lines present the dependence of 
[Oin] A5007/H/3 (upper panel), [Oin] A4363/[Oni] A5007(middle 
panel), and H/3 EW (lower panel) on n^. at given values of U, 
as calculated by CLOUDY. The square boxes represents objects 
with detections, and the triangles objects with an upper limit 
on the [Oiii] A4363 EW, i.e. an upper limit on Ue- In the lower 
panel the calculated H/3 EW is for CF = 1, and for each object 
we therefore plotted the observed H/3 EW divided by the inferred 
CF. 



line strengths and profiles (subject to the condition that 
broad emission lines are present). 

Optical observations of the H/3 region of the 87 objects 
are described in BG92, and were kindly provided by T. Boro- 
son (private communication). These spectra do not extend 
down to [O III] A4363 in all objects, and we therefore sup- 
plemented the EG spectra by optical spectra from Shang 
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Figure 3. The theoretical and measured relation between 
[Oiii] A5007/H/J and [Ooi] A4363/[Oiii] A5007. The solid lines 
represent curves of constant Ue and constant U, as calculated by 
CLOUDY. The object symbols are as in Fig. 2. 



et al. (2003), kindly provided by Z. Shang (private commu- 
nication) which also provide a higher S/N in some cases'^. 
The final sample contains 78 objects with spectra which con- 
tain the [O III] A4363 region. The spectra were corrected for 
reddening, redshift, and possible slit losses, as described in 
Baskin & Laor (2005, hereafter BL05, Section 2.2 there). 



3.2 The Fell subtraction 

The [O III] A4363 line is generally very weak, and its mea- 
surement requires a careful subtraction of the Fe ll emission 
multiplets. We used the 1 Zw 1 Feii template, kindly pro- 
vided by T. Boroson, to subtract the Fell lines from the 
spectra (see Section 2.2 in BL05 for further details concern- 
ing the subtraction procedure). 

The Fell template was constructed by subtracting all 
the none Fe ll emission from the observed spectrum of 1 Zw 1. 
However, this also requires estimates of the relative contri- 
butions of Fe II and [O ill] A4363 near 4363A, which may 
not be accurate, and could lead to biased estimates of the 
[O III] A4363 strength in the Fe ii subtracted spectra. To 
check the accuracy of the Feii subtraction, we also used 
the more theoretically motivated Fen template of 1 Zw 1, 
independently constructed by Vcron- Getty, Joly & Veron 
(2004), kindly provided by M. Veron-Cetty. The BG92 tem- 
plate is rather smooth and featureless around 4363A, and 
the Veron-Getty et al. template has a sharp feature redward 
of 4363A, but much weaker emission at 4363A, compared 
to BG92. The Veron-Getty et al. compilation of expected 
allowed and forbidden Fell emission lines indicates no sig- 
nificant feature near 4363 A. There was almost no effect on 
the [Oiii] A4363 profile in 1 Zw 1 when using either the 
Veron-Getty et al. template or the BG92 template to sub- 
tract the Fe ii emission. We therefore conclude that improper 

^ The spectra obtained from Shang et al. (2003) are of 
PG 0953-1-414, PG 1216-1-069, PG 1322+659, PG 1427-1-480, 
PG 1512-f370 and PG 1543-1-489. 
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Fe II subtraction is unlikely to significantly bias our measure- 
ments of the [O III] A4363 line. In our analysis we used the 
empirically derived BG92 Fe ii template. 



3.3 The measurement procedure 

As noted above, there is likely to exist a distribution of Ue, 
U, and CF values for the [Oiii] emitting gas in the NLR. 
To ensure as much as possible that the measured H/3 and 
[O III] A4363 emission originates in the same gas which pro- 
duces the [O III] A5007 emission, we use the [O ill] A5007 
profile as a template for measuring the H/3 and [O ill] A4363 
lines. This ensures that all three lines originate in gas with 
the same line of sight velocity distribution, although it does 
not necessarily guarantee they all originate in NLR gas 
clouds with the same distribution of emission properties. 

In each object we subtracted the ma:ximum allowed 
scaled [O ill] A5007 profile which does not produce a dip in 
the broad H5 line and in the H7-|-[0iii] A4363 blend, as de- 
termined by eye inspection. The subtraction was performed 
as follows. Both the [O ill] A5007 line and either [O ill] A4363 
or H/3 were transformed from wavelength scale to velocity 
scale, using z measured from the [O ill] A5007 peak (T. Boro- 
son, private communication, listed in BL05). Then, the max- 
imal c in the term /linc(t^) — cx /5007(i') was searched, where 
'line' is either H/3 or [O ill] A4363, such that the subtrac- 
tion of the scaled [O iii] A5007 profile will not produce a dip 
in the measured line. The measured narrow line component 
equals the scaled [Olll] A5007 line, /linG(i') = c x fsoo7{v)- 
To calculate the line EW the line profile func{v) was trans- 
formed back to wavelength dependence /line(A), integrated 
over A, and divided by the continuum flux density at 4861 A 
(the same wavelength used by BG92). The flux density at 
4861 A was measured by fitting a local power-law contin- 
uum to each spectrum between ^ 4600 A and ^ 5100 A. 
We measured the line EW directly, rather than infer it from 
the [O III] A5007 EW measured by BG92 and the value of c, 
since the [O iii] A5007 profile deduced here may not be iden- 
tical to the one assumed by BG92 due to possibly different 
definitions of the continuum underlying [O ill] A5007. 

To determine whether the measured line flux is signif- 
icant or not, we need to evaluate the standard deviation, 
<T, of the flux measurement. We flrst measured the RMS of 
the flux density between 4500A and 4600A^, multiplied the 
RMS by the wavelength width of a pixel (to convert flux 
density to flux), and then multiplied it by the square- root 
of the number of pixels in the full width at zero intensity 
of the [O III] A5007 profile, to get a. A measured line flux 
was accepted as significant if it is > 3a. Following BG92, we 
also adopted an EW of 0.5 A as the lower detection limit 
(accounting for the minimum likely level of systematic er- 
rors). After exclusion of objects with an unmeasurable H/3 
narrow component, and an unmeasurable [O ill] A4363, the 
remaining sample contained 40 objects where all three lines 
are detected. 

Figure 4 shows four representative cases of the 
[Oiii] A5007 profile subtraction from the H7+P111] A4363 
blend. The two left panels present objects where the 



^ Seven spectra have a feature in this range, and thus a slightly 
different range was chosen for them. 
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Figure 4. The H7 and [O in] A4363 blend before (top con- 
tinuous line) and after (dashed line) substraction of the scaled 
[Oiii] A5007 profile, and the [Oiii] A5007 profile (scaled by an 
arbitrary factor for presentation; bottom continuous line). The 
name of the object is indicated in each panel. Upper-left, the 
significance of the measurement is 18.6<T. Lower-left, the signifi- 
cance is 3. 60". The two right panels demonstrate a broad residual 
[O III] A4363 component that often remains (see text) . 



[O III] A4363 emission is rather well fit using the [O in] A5007 
profile. The upper left panel is a clear 18.6cr detection, and 
the lower left one a marginally acceptable detection at 3.6a. 
The two right panels in Fig. 2 demonstrate cases where a 
broad residual [Oiii] A4363 component remains after sub- 
traction of the scaled [Oiii] A5007 profile. Such a broad 
residual component often remains, and it most likely indi- 
cates an inner higher density extension of the NLR, where 
the density is higher than the [O ill] A5007 critical density 
but it is lower than the critical density of [Oiii] A4363. Such 
an inner, higher velocity, extension of the NLR would con- 
tribute to the [O III] A4363 line and not to the [O ill] A5007 
line (Nagao et al. 2001). This demonstrates that making in- 
ferences based on the total [Oiii] A5007 and [Oiii] A4363 
line fluxes, rather than matched proflle fluxes, could lead to 
inaccurate results. 



3.4 The inferred NLR parameters 

Table 1 presents the narrow "[Olll] A5007-like" H/3 and 
[O III] A4363 EW, measured as described above. We also list 
in Table 1 the [O ill] A5007 EW from BG92 (the Shang et al. 
spectra yielded consistent values). Column 1 lists the name 
of the object, columns 2 and 3 list the measured H/3 EW in 
units of A and its significance in units of a. Column 4 lists 
the [O III] A5007 EW in units of A, as reported by BG92. 
Columns 5 and 6 list the measured [O iii] A4363 EW and 
its significance. Column 7 lists log vL,, at 4861 A in units of 
10'*'' erg s~^, which allows a quick conversion of line EW to 
line luminosities. Columns 8, 9 and 10 list the inferred NLR 
parameters, logt/", logrie (in units of cm^'^), and CF (in per 
cent) for the single-zone model. Finally, columns 11, 12 and 
13 fist the CFin (in per cent), logUout and CFout (in per 
cent) for the two-zones model. 

To infer U, He and CF we calculated theoreti- 
cal values of the H/3 flux, and [Oiii] A5007/H/3 and 
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Table 1. The measured emission line parameters and the inferred NLR parameters. 



Object 


H/3 




[Oiii] 


[Oiii] 


A4363 




Sin 


gle zone 




Two zones 






EW 


S/N 


A5007 


EW 


S/N 




U 


n 


CF 






CFout 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 




1.96 


20.9 


26 


1.24 


16.5 


1.579 


—2.2 


5.64 


3.6 


4.8 


— 1.50 


4.7 




5.34 


22.9 


22 


1.04 


6.8 


—0.267 


—3.3 


5.79 


6.8 


4.8 


—3.18 


6.1 


0007-4-1 Ofi 


4.19 


16 


42 


0.84 


4.4 


0.407 


—2.9 


5.24 


5.7 


2.8 


—2.45 


5.8 




3.10 


34.9 


29 


1.54 


21.2 


0.804 


—2.6 


5.75 


4.6 


6.9 


—2.23 


3.9 




15.55 


18.8 


97 


6.07 


10.9 


—0.509 


—3.0 


5.88 


20.5 


28.9 


—2.85 


16.4 


on "^0-1-1 94 


Q7 


41 4 


22 


4 19 




n "^1 9 

U.O J.Z 





6.40 


9.2 


21 4 


3 22 


^ Pi 
0.0 




3.90 


29.9 


30 


1.22 


10.3 


0.808 


— ^ 


5.67 


5.2 


5.4 


—2.68 


4.6 




3.58 


8.8 


13 


1.40 


4.4 


0.356 


—3.1 


6.15 


4.7 


7.0 


—3.20 


3.5 


0091 -l-'^9'^ 


13.00 


33.3 


83 


5.46 


18.6 


—0.659 


—3.0 


5.91 


17.2 


26.0 


—2.80 


13.5 


uyzon^ ± zy 


7.04 


24.7 


31 


1.56 


8.2 


—0 579 


—3.2 


5.81 


9.1 


7.3 


—3.15 


8.1 


u y o'-in^ u J- 


7.29 


22.6 


51 


2.58 


11.7 


—0.603 


—3.0 


5.78 


9.6 


11.9 


—2.75 


8.2 


VjyOO^'i -L'i 


1.87 


10.0 


18 


1.16 


7.2 


1.033 


—2.6 


5.85 


2.8 


5.3 


—2.01 


2.3 


1 001 -UO'^zL 


2.65 


5.4 


7 


1.33 


3.2 


0.480 


—3.1 


6.41 


3.4 


6.8 


—3.34 


2.3 


1 048-1- '^42 


2.72 


5.9 


28 


1.90 


5.5 


1.235 


—2.4 


5.85 


4.4 


8.4 


— 1.71 


4.4 


1 1 00-1-779 


3.90 


15.8 


41 


0.98 


4.4 


1.276 


—2.8 


5.34 


5.4 


3.5 


—2.32 


5.6 


111 Q-Ll 90 


2.66 


23.0 


19 


0.91 


9.7 


—0.140 


—3.0 


5.76 


3.5 


4.2 


—2.73 


3.0 


1 1 9R— 041 

X i.£i\J — utx 


3.47 


12.4 


19 


0.60 


3.2 


0.021 


—3.2 


5.57 


4.5 


2.6 


—3.05 


4.1 


1149—110 


4.64 


16.7 


33 


2.07 


9.5 


—0.392 


—2.9 


5.87 


6.3 


9.8 


—2.66 


4.9 


1151+117 


1.33 


17.1 


w 


0.66 


10.6 


0.456 


—2.8 


5.84 


1.9 


3.1 


—2.44 


1.5 


1 2(12-1-2X1 

_L ij X 


5.64 


37.3 


36 


1.30 


9.8 


0.386 


—3.1 


5.63 


7.4 


5.7 


—2.93 


6.6 


1 99Q_L904 


1.94 


10.9 


19 


0.79 


6.5 


0.067 


—2.7 


5.64 


2.8 


3.4 


—2.26 


2.6 


_L irf rt rt 1^ V/ ^ V/ 


8.21 


28.1 


17 


0.75 


3.5 


—0.454 


—3.5 


5.78 


10.2 


3.5 


—3.51 


9.6 


1 Q09— 1 09 


2.16 


12.8 


9 


0.88 


5.5 


1.510 


—3.1 


6.11 


2.8 


4 4 


—3.13 


2 1 




3.18 


14.7 


32 


1.80 


11.8 


0.716 


—2.6 


5.78 


4.7 


8.0 


—2.01 


4.2 




2.55 


18.0 


19 


0.89 


7.1 


0.756 


—3.0 


5.75 


3.4 


4.1 


—2.68 


2.9 


1 mo— 1 v\R 

xoxu — xuo 


8.24 


39.2 


77 


5.41 


30.6 


—0.637 


—2.5 


5.88 


12.5 


24.8 


—2.01 


10.0 


1 ^41 4-9'^« 


1.90 


9.8 


14 


0.54 


4.1 


—0.053 


—3.0 


5.65 


2.5 


2.4 


—2.73 


2.2 


1 ^^1 4-9^fi 


2.46 


11.2 


12 


0.57 


3.0 


—0.204 


— ^ 9 


5.78 


3.2 


2.6 


—3.09 


2.8 


1351-1-640 


3.70 


20.6 


31 


1.77 


11.3 


0.515 


-2.8 


5.81 


5.2 


8.2 


-2.44 


4.1 


1352+183 


1.57 


4.8 


10 


1.33 


3.6 


0.432 


-2.6 


6.19 


2.3 


6.6 


-2.34 


1.2 


1411+442 


2.25 


9.9 


15 


1.01 


6.5 


0.279 


-2.9 


5.91 


3.1 


4.8 


-2.73 


2.3 


1425+267 


3.61 


20.0 


36 


0.84 


4.9 


1.301 


-2.8 


5.33 


5.0 


3.0 


-2.42 


5.0 


1501+106 


5.58 


35.6 


64 


1.20 


10.9 


0.016 


-2.7 


5.15 


8.0 


3.6 


-2.22 


8.8 


1512+370 


5.50 


20.9 


57 


1.14 


5.2 


0.958 


-2.8 


5.22 


7.7 


3.8 


-2.39 


7.9 


1519+226 


5.87 


8.6 


4 


2.45 


4.8 


0.265 


-3.2 


7.04 


7.5 


12.8 


-3.96 


5.4 


1534+580 


5.99 


48.3 


79 


2.33 


25.3 


-0.542 


-2.5 


5.42 


9.0 


7.8 


-1.75 


12.5 


1535+547 


4.53 


22.0 


16 


0.52 


4.3 


-0.286 


-3.4 


5.61 


5.7 


2.3 


-3.27 


5.4 


1545+210 


2.61 


11.2 


33 


1.20 


6.5 


1.093 


-2.5 


5.54 


4.0 


4.4 


-1.72 


5.3 


1612+261 


14.28 


74.3 


157 


3.11 


23.1 


0.332 


-2.8 


5.22 


19.9 


9.9 


-2.28 


21.6 


1704+608 


2.70 


15.2 


27 


0.52 


3.4 


1.410 


-2.9 


5.23 


3.7 


1.8 


-2.46 


3.7 



[Oiii] A4363/[Oiii] A5007 ratios using the photoionization 
code CLOUDY (v. 94.00). The calculations wore made for 
a grid of models with —4 < log U < — 1 (in steps of 0.5), and 
3 < log Tie < 8 (in steps of 0.1). The H/3 flux was converted 
to EW as follows. The choice of U and Ue corresponds to a 
given ionizing flux at the face of the photoionized slab, which 
is specified by the output of CLOUDY. We multiplied the 
H ionizing flux by 0.4 to get an approximate value for vfi, 
at A4861 (deduced from the spectral shape in fig. 7 of LD), 
and divided it by 4861 to get the continuum fx at 4861 A. 
The H/3 fine flux from CLOUDY was then divided by this 
fx to get the predicted line EW for CF — 1 . This procedure 
assumes a uniform spectral energy distribution (SED) in all 
objects, which may not correspond well with the range of ob- 
served SEDs in the PG sample quasars. However, as shown 
in Section 4.1.3, the likely range of SED shapes has only a 
minor effect on the deduced parameters. 

For the single zone model we use the following algorithm 



to infer U, Ue and CF from the H/3 EW, [O iii] A5007/H/3 and 
[Oiii] A4363/[Oiii] A5007 ratios. 

(i) Assume U, starting at log U = —4. 

(ii) Calculate rie from the observed 
[Oiii] A4363/[Oiii] A5007 ratio for the assumed U 
(Fig. 2, middle panel). The calculation is made using a 
linear interpolation scheme for log [O iii] A4363/[0 iii] A5007 
versus log rie, when the observed [Oiii] A4363/[Oiii] A5007 
value is between grid points. 

(iii) Calculate the theoretical [Olll] A5007/H/3 ratio for 
the assumed U, and the Ue calculated above (using a linear 
interpolation scheme as above). 

(iv) Calculate the theoretical [O ill] A5007 EW from the 
theoretical [O iii] A5007/H/3 ratio found above, and the ob- 
served H/3 EW. 

(v) Repeat steps (i)-(iv) changing \ogU from —4 to — 1 
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in steps of 0.1 (using linear interpolations between U grid 
points). 

(vi) Adopt as the solution the values of U and rie which 
yield the smallest difference between the observed and the- 
oretical [O III] A5007 EW, as calculated in step (iv) . 

(vii) Calculate the theoretical H/3 EW for the U and rie 
of the solution (Fig. 2, lower panel). Obtain CF by dividing 
the observed H/? EW by the theoretical one. 

This procedure is effectively identical to the one described in 
Section 2. The differences between the observed and calcu- 
lated [O III] A5007 EW are typically well below 10 per cent, 
which is generally within the [O ill] A5007 EW measurement 
errors. The resulting best fit U, Ue and CF for each object 
are listed in Table 1. 

As discussed in Section 2, in the two zone model 
we assume Ue = W'^ cm^'^ for the outer zone, and 
He = lO'^ cm~^ and U = 0.1 for the inner zone, leav- 
ing CFin, CFout and log?7out as the remaining three pa- 
rameters to be sot by the observations. The assumed 
and U for the inner zone imply log [Olll] A5007/H/3=0.2, 
log [Olll] A4363/[Olll] A5007=-0.04 (see Fig. 3) and a 
calculated H/3 EW of 12.21 A for CFin=l. We now iter- 
ate over the possible range of log U values for the outer 
zone. For each U value for the outer zone we calculate 
[Oiii] A5007/H/3 [O iii] A4363/[0 iii] A5007, and the H/3 EW 
for CFout = 1. Wc now write two linear equations for the ob- 
served [O III] A5007 and [O iii] A4363 EWs in a given ob- 
ject, as a function of the above line ratios, the calculated 
H/3 EWs, and the two unknowns, CFin and CFout, which 
we solve for. The solution produces the correct observed 
[Out] A5007/H/3 [Oiii] A4363/[Oiii] A5007 ratios by con- 
struction, and it also allows us to calculate the total ex- 
pected H/3 EW. The accuracy of the solution is measured 
by the difference between the observed and calculated total 
H/3 EW. The U value which minimizes this difference (typ- 
ically less than 10 per cent) is taken as the best fit solution 
for the given object. 



4 RESULTS &: DISCUSSION 

4.1 The distribution of rie and U values 

4.. 1.1 The Single Zone Approximation 

As shown in Table 1 and Fig. 2, the values of rie of all 40 ob- 
jects with detections (excluding one outlier, PC 1519-1-226) 
lie within a relatively small range (factor of ~ 5) around 
the critical density of [Olll] A5007 (i.e. within log rie = 
5.85±0.7). This is expected if there is a range of densities in 
the NLR, since the line emissivity per unit mass peaks close 
to the critical density, and thus the [O iii] A5007 emissivity 
weighted mean density is expected to be at Ue ~ ricrit, un- 
less the distribution of CF with density is heavily weighted 
towards very high or very low densities. 

We also find that all objects (except one) lie at log U < 
—2.2. A similar selection effect may be present here as well. 
As Fig. 2 (lower panel) shows, the H/3 line emissivity per unit 
solid angle (= EW/CF) drops at logL*' > -2.5 (due to dust 
suppression), and thus if there is a distribution of U values 
in the NLR with comparable values of CF, most of the H/3 
emission would originate in gas with logZ7 < —2. However, 



at log U < —3.5 the O ill emission is suppressed because 
Oxygen becomes less ionized than O ill. Thus, the combined 
effects of dust suppression and ion abundance leads to a 
relatively small range of U (log U = —3 ± 0.5) where both 
the [O III] lines and H/3 are emitted efficiently. 

In the above analysis we considered only objects where 
all the three narrow lines are detected. Ignoring upper limits 
may result in spurious trends or correlations. To explore 
that, we repeated the analysis described in Section 3 for all 
the objects where [O ill] A4363 was not detected, assuming 
the line EW equals the upper limit value (the larger of the 
3cr detection limit and 0.5 A). These objects are also plotted 
in Figs. 2 and 3, and their distribution nearly overlaps the 
distributions of the detected objects. We also repeated the 
analysis arbitrarily assuming the [O ill] A4363 EW is a factor 
of 10 lower than the upper limit. This leads to a tail of 
objects extending to log rie = 3, but the distribution of U 
values is not much affected. The extension to low rie occurs 
because rie is largely set by [O ill] A4363/[0 ill] A5007, while 
U is unchanged because it is largely set by [O ill] A5007/H/3, 
and it is nearly independent of the [Olll] A4363 EW when 
[Oiii] A4363/H/3 is low (Fig. 3). 

4- 1.2 The Two Zone Approximation 

As described above, the two zone approximation we use as- 
sumes fixed values for rie in the two zones, and a fixed U for 
the inner zone. For [Tout we find —4 < log C/out < —1.5, with 
most objects in the range of —3.5 to —2. This distribution 
is offset to somewhat higher U values (by ^ 0.2) compared 
to the single zone approximation (Fig. 5). 

4.1.3 Modified SED 

Far UV observations of high z quasars (Zheng et al. 1997), 
together with soft X-ray observations of low z quasars (Laor 
et al. 1997), suggest a significantly softer ionizing SED than 
the one suggested by MF, which we used in the above anal- 
ysis (but see Scott et al. 2004). To explore the possible role 
of the SED on the deduced NLR parameters we repeated 
the analysis in Section 3 with a modified softer SED, where 
we added a break at 1 Rydberg to a power-law slope of —1.5 
extending to 70 Rydberg {r^ 1 keV), as suggested by the far 
UV to soft X-ray composite of Laor et al. (1997). The rest 
of the SED remains unchanged. 

The resulting softer SED values, nl, and If, are very 
similar to those obtained with the harder MF SED. Apart 
from two outliers, where log U changed from ^ —2.3 to 
~ —1, due to the degeneracy of the U solution at high U 
(Fig. 3), for the other objects we find a very tight relation 
between U and U" . Specifically, we find a best fit relation 
log [/" = 0.9 log U - 0.24, with an RMS scatter of 0.05, and 
a Spearman rank order correlation coefficient rs = 0.989. 
This relation implies an average systematic ofl^set of log U" — 
log[/ = -0.04 to 0.11, for logU = -2 to -3.5, which is 
well within the likely systematic and statistical errors in our 
solution for U. Similarly, for nl we find log = 1.02 log rie — 
0.17, an RMS scatter of 0.02, rs = 0.996 (excluding the 
two outliers), and again a negligibly small offset of 0.07 to 
0.03 for log rie = 5 to 7. For the soft SED covering factor 
CF" we get CF" = l.OOCF-l-0.24, with an RMS scatter of 
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only 0.01 and rs = 0.999. Thus, CF' is essentially perfectly 
correlated with CF, but it is offset to higher values by a 
factor of 1.74(=10"'^''), as expected since the softer SED 
requires a larger CF to produce a given H/3 EW. Note that 
the exact absolute value of the best fit parameters, rie, U, 
and CF, for each object are not crucial for our study of the 
modulation of the [Oiii] A5007 strength, and what mostly 
matters is the range of values for these parameters. 



4-i-4 Comparison with Earlier Studies 

Binette et al. (1996, and references therein) discuss evidence 
for the presence of significant contribution to the NLR emis- 
sion from matter bounded gas. This component produces 
mostly high ionization lines, and it is invoked to explain the 
large relative strengths of the Hell A4686 and CiV A1549 
lines, seen in the NLR emission of Seyfert 2 galaxies (e.g. 
Ferland & Osterbrock 1986), as well as the strength of 
other lines not measured in our sample. The single and two 
zone models we use do not include such a matter bounded 
gas component. However, we suspect that such a compo- 
nent may not be required by our data, and the anomalous 
strength of Hell A4686 may be interpreted differently, as 
further discussed below. 

Binette et al. (1996) used a mean Civ A1549/H/3 flux 
ratio of 12 ±8 from Ferland & Osterbrock (1986). However, 
using the narrow Civ A 1549 EW measurements available 
for our sample (BL05) , wc deduce a significantly lower mean 
value of 4.7 ± 6, which is likely biased to large values as it 
does not exclude low significance level narrow CiV A1549 
measurements. The weakness of the NLR CiV A1549 com- 
ponent in broad line AGN was already noted by Wills et al. 
(1993) in Hubble Space Telescope {HST) spectra of seven 
radio loud AGN*. The lower mean CiV A1549/H/3 value is 
consistent with standard photoionization model predictions 
(e.g. Wills et al. 1993; Groves et al. 2004a), and thus there 
is no evidence from C iv A1549 for a matter bounded NLR 
component in our objects. 

The mean He ii A4686/H/3 flux ration in Seyfert 2 galax- 
ies is 0.29 ± 0.1 (Ferland & Osterbrock, 1986). We measure 
a similar mean value of 0.30 ± 0.12 in 43 of our objects 
where a narrow "[O ill] A5007-like" Hen A4686 is detectable. 
As shown in the extensive set of calculations by Groves 
et al. (2004a; 2004b), dusty photoionized gas models with 
log [/ ~ —2 can produce Hell A4686/H/3~ 0.3 in ionization 
bounded gas. The models of Binette et al. (1996) do not 
include dust, and these produce lower Hell A4686/H/3 ra- 
tios, which lead Binette et al. to conclude that an additional 
population of matter bounded gas clouds must be present at 
the NLR. Dusty photoionized gas with log U ^ —2 produces 
a higher Hell A4686/H/3 ratio since the H/3 line originates 
from deeper layers compared to Hell A4686, and is thus more 
strongly suppresed by dust absorption. Most of the line emis- 
sion in dusty photoionized gas originates within a surface 
layer with a column of < 10^^ cm^'^, where the dust opacity 
is < 1. The line emission from this surface layer mimics to 



* The apparently discrepant C iv A1549 NLR emission in Type 
II/Type I AGN is inconsistent with the inclination unification 
schemes for AGN, and this needs to be addressed in future studies. 



some extent the line emission from matter bounded, dust 
free photoionized g invoked by Binette et al. (1996). 

Dopita et al. (2002) calculated photoionization mod- 
els for isobaric dusty gas, and proposed that the ionization 
parameter in such gas tends to a constant value, indepen- 
dent of the surface gas density, once radiation pressure dom- 
inates gas presure. Dopita et al. further suggested this ef- 
fect explains the small dispersion in the NLR line ratios 
(log [Oiii] A5007/H/3~0.9-1.2 in their sample), and the im- 
plied nearly constant \ogU ^ —2 in the NLR. However, 
here we flnd a much larger range (log [Oiii] A5007/H/9~0- 
1.1, Fig. 3), which implies a correspondingly larger range in 
U. This argues against the effectiveness of the mechanism 
proposed by Dopita et al. to produce a fine tuned U in the 
NLR. We note in passing that the dusty isochoric CLOUDY 
model used here produces nearly identical line ratios to 
those reported by Groves et al. (2004a; 2004b) for isobaric 
dusty gas with significant radiation pressure. For example. 
Groves et al. find log [Olll] A4363/[Olll] A5007= -1.65 for 
rie = 10^ cm~^, U = 0.1, and solar abundance (Fig. 13 in 
Groves et al. 2004b), while we find —1.75 for the same pa- 
rameters (Fig. 3). At U = 10~*, where radiation presure 
effects on dust are negligible, their model gives a ratio of 
—2.3, and CLOUDY gives an essentially identical value of 
-2.36. 

4.2 Which parameter modulates most strongly 
the [O III] A5007 EW? 

Figure 5 (upper panels) presents the dependence of the 
[O III] A5007 EW on U, Ue , and CF, for the single zone model. 
The [O III] A5007 EW ranges over a factor of 40 (4A to 157A) 
in our sample of 40 AGN with complete detections (com- 
pared with a range of > 300 in the complete sample of 87 
objects). The [O III] A5007 EW is most strongly correlated 
with the CF. The Spearman rank order correlation cocffl- 
cient, rs, for the 40 detected objects is 0.699, which has a 
null probability of Pr=5.2 x 10~^ (note that the distribution 
of the [O 111] A4363 non detections is consistent with the dis- 
tribution of the detections) . The CF covers a range of about 
10 (1.9 to 20.5 per cent. Table 1), and since the EW of aU 
lines is linear with CF, the CF modulates the EW by a fac- 
tor of 10, making it the dominant factor in modulating the 
[Oiii] A5007 EW. 

The [Oiii] A5007 EW is also significantly correlated 
with Tie {rs = -0.473, Pr = 2 x 10"^), though the 
[O III] A4363 non detections may affect the strength of this 
correlation if the true [O in] A4363 EW are well below the 
upper limits. The range of rie ^ 10^ — 10®'^ cm^"'' amounts to 
a modulation of the [O ill] A5007 EW by an additional factor 
of 3 (e.g. the log U — —2.5 curve in Fig. 2, upper panel). Fi- 
nally, the [O III] A5007 EW shows only a marginal correlation 
with U {rs = -0.384, Pr = 1.5 x 10"^, not affected by upper 
limits). Although the range of C/ ~ lO-^ '^-lO"^'^ modulates 
the [O III] A5007 EW by a factor of ~ 4, its weak correlation 
with [O III] A5007 EW indicates it mostly contributes to the 
scatter in the previous two relations, and does not have a 
large systematic effect on the range of [O III] A5007 EWs. 

The lower panel of Fig. 5 presents the dependence of 
the [Oiii] A5007 EW/CF on U,ne and CF. The [Oiii] A5007 
EW/CF represents the [O in] A5007 emissivity per unit solid 
angle of the NLR gas. This parameter ranges over a fac- 
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tor of 5, which represents the amount of [O ill] A5007 EW 
modulation by the combined effects of rie and U. The NLR 
[O in] A5007 emissivity per unit sohd angle is strongly corre- 
lated with U {rs = -0.800, Pr = 5.9 x 10"^°) and somewhat 
less strongly with Ue {rs ~ —0.708, Pr = 3.2 x 10^^, again 
weakened by non detections). This emissivity is completely 
unrelated with CF (rs = -0.034, Pr = 0.83). We finally note 
that He, U, and CF, are not significantly correlated with each 
other. 

As expected (Section 4.1.3), essentially identical re- 
sults were obtained with the modified soft SED. Somewhat 
stronger correlations were obtained in the two zone approxi- 
mation. Specifically, the [O ill] A5007 EW correlation with U 
rose to 0.454 and with CFout to 0.778 (both parameters for 
the outer zone), while the correlations for the [Oiii] A5007 
emissivity (Fig. 5, lower panel) remained essentially un- 
changed. 

4.3 Relation with EVl parameters 

The [O III] A5007 line is one of the main components of the 
BG92 EVl set of correlations. These correlations may be 
driven by some fundamental parameters, such as L/I/Bdd 
(BG92) . Below we briefiy explore whether the NLR physical 
parameters ne,U, and CF are correlated with some of the 
non- [Oiii] A5007 EVl components (as listed in BG92). 

The CF is significantly correlated with the absolute V 
band magnitude Mv {rs = 0.512, Pr = 7.6 x 10""), such 
that the CF decreases with increasing luminosity. This ac- 
counts for the tendency of the [O ill] A5007 EW to decrease 
with increasing luminosity in the BG92 sample (note that 
Shcmmcr ct al. 2004 do not find this efi^cct in luminous liigli- 
z quasars) . The only other parameter significantly correlated 
(i.e. Pr < 10"^) with CF is the EW of the broad He II A4686 
line {rs = 0.452, Pr = 3.4 x 10~*). This correlation is not 
simple to interpret, and it may be induced by the strong 
inverse correlation of the broad He II A4686 EW with lumi- 
nosity, and the above inverse correlation of luminosity and 
CF. Interestingly, there is no significant correlation of CF 
and L/Ljsdd (taken from BL05; rs = -0.133, Pr = 0.41). 
Such a relation was suggested by BG92 as the physical mech- 
anism driving the set of EVl correlations (see also Kraemer 
et al. 2004). 

The rie is significantly correlated with two surprising pa- 
rameters, the soft X-ray slope (rs = —0.633, Pr = 1 x 10"''), 
and the compact to total radio flux ratio from Kellermann 
et al. (1989) (rs = 0.560, Pr = 1.7 x 10""). A marginally 
significant correlation exists with the Fe II/H/3 flux ratio 
(rs = 0.415, Pr = 7.8 x 10"^), which may point at a rela- 
tion between the NLR density and the BLR metalicity. 

The U is significantly correlated with the continuum 
luminosity at 3000 A, and with the broad Rjd line FWHM, 
leading to a somewhat stronger, but hard to interpret corre- 
lation (rs = 0.552, Pr = 2.3 x 10"") with the estimated black 
hole mass (Mbh, sec BL05). A stronger correlation exists 
with the Fe II/H/3 fiux ratio (rs = -0.608, Pr = 3.1 x 10"^). 
Again, this may point at a relation between the NLR ion- 
ization level and the BLR metalicity. 

The two zone analysis yields similar correlations for 
CFout with Mv {rs = 0.456, Pr = 3.1 x 10"^), and with 
the broad He II A4686 EW (rs = 0.460, Pr = 2.8 x 10"*). 
But, there are no significant correlations with CFin. The cor- 



relations for f/out arc somewhat stronger than those for the 
single zone U (with A/bh, rs = 0.590, Pr = 6.1 x 10"^; with 
Fe II/H/3 rs = 0.654, Pr = 4.8 x lO"''). 

4.4 The size of the NLR 

The distance of the [O lll] A5007 emitting region from the 
central ionizing continuum source, -Rnlr, can be determined 
from the measured rie and U for each object, and the esti- 
mated ionizing luminosity Lion, as follows. By definition U = 
n-i/ue, where the ionizing photon density is given by n-y = 
i'ion/47ri?NLRc(fei^), where (hv) is the mean ionizing pho- 
ton energy (3.03 Rydberg, for the CLOUDY MF ionizing 
continuum). These relations give /^nlr = KL\/^^{Une)~^^^ , 
where K = (47rc(fti/))"i/2. 

Figure 6 presents the luminosity dependence of -Rnlr 
for the single zone and for the two zone approximations. 
The correlation for the single zone is strong (rs = 0.836, 
Pr = 2 X 10^^^), with a best fitting relation 

-Rnlr = 40-L44"^ pc, 

where L44 = !vL^(4861)/10""erg s"\ For the two zone ap- 
proximation we get a somewhat weaker but still signifi- 
cant correlation of the size of the outer zone (rs = 0.648, 
Pr = 6 X 10"®), with a best fitting relation 

7-)Out Ten r 0.34 _ 

rtNLR = 75O-L44 pc. 

The assumption of fixed ne and U values in the inner zone 
implies -Rnlr = ^44^ PC by construction. 

The single zone -Rnlr versus L relation implies sizes 
which are much smaller than deduced through imaging of 
spatially extended [Oiii] A5007 emission in Type I AGN 
with the HST (Bennert et al. 2002; Schmitt et al. 2003). 
For example, in four overlapping PC quasars (0026-1-129; 
0053-1-251; 0953-^414; 1307-^085) Bennert et al. find i?NLR = 
2.2-5.9 kpc, while we get .Rnlr = 84-116 pc, i.e. about ~ 40 
times smaller sizes. In the two zone approximation we get 
^NLR =1.3-1.7 kpc, i.e. within a factor of 2-3 of the Bcrmcrt 
et al. values. However, the Bennert et al. HST images ap- 
pear to be measuring the kpc-scale 'extended NLR' (ENLR, 
Unger et al. 1987), which according to the surface bright- 
ness values in Bennert ct al. (Table 2 there), may include a 
significant component of unresolved [O ill] A5007 emission. 
A more carefuU analysis of these images is required to deter- 
mine the fraction of unresolved [O iii] A5007 emission, and 
the implied constraints on the [O ill] A5007 radial emissivity 
distribution. Note also that the extrapolation of the Ben- 
nert et al. relation to higher L implies -Rnlr values which 
are ruled out by ground based observations of high z quasars 
(Netzer et al. 2004), also pointing at a more compact .Rnlr. 
Similarly, a compact NLR (< 100 pc), consistent with our 
estimates, was deduced by Kraemer et al. (1998), based on 
detailed photoionization modeling of the narrow line emis- 
sion in NGC 5548. 

The slope of 0.45 for the single zone model implies that 
Z7ne, i.e. is roughly constant. The scatter in the -Rnlr 
versus L relation is only 0.2 (RMS in log -Rnlr), which is 
significantly smaller than expected based on the measured 
range of values for U and for Ue ^ ■ The small spread in the 

® Phenomenologically, this results from a weak anti-correlation 
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Figure 5. The dependence of [Oiii] A5007 EW (three upper panels) and [Oiii] A5007 EW/covering factor (three lower panels) on the 
ionization parameter, density and covering factor. The Spearman rank-order correlation coefficient (rg) and the null probability (Pr) are 
indicated in each panel. 



Anlr versus L relation, and a slope close to 0.5, indicate a 
small spread of n-, at the [O ill] A5007 emitting region. It is 
not clear why such a small range in n-y should be present 
The two zone model shows a somewhat larger scatter (0.3 
RMS in log i?NLR)i and a flatter slope (0.34), which indicates 
a non-uniform n~f. 

The tight i?NLR versus L relation may be partly induced 
by the rejection of objects with non-detected [O ill] A4363. 
In these objects, the arbitrary assumption that the true 
[Oiii] A4363 EW is 0.1 times the upper limit value leads 
to i?NLR values which are an order of magnitude larger than 
deduced for the rest of the objects. Thus, higher quality op- 
tical spectra, with complete detections of [O iii] A4363, are 
required in order to measure more accurately the amount of 
scatter in the 7?nlr versus L relation. 

Figure 7 shows that there is also a significant in- 
verse correlation between 7?nlr and rie (rs ~ —0.596, 
Pr = 4.9 X 10"^). The slope of the relation is -0.56. This re- 
lation is consistent with earlier findings of a negative radial 



density gradient in the NLR. In this case, inclusion of upper 
limits on [O ill] A4363 as detections, at a tenth of their value, 
extends the correlation to large 7?nlr and lower rie, indicat- 
ing that this correlation cannot be induced by ignoring the 
non detections. 

Inspection of the correlations of J?nlr with the EVl 
parameters yields that the broad H/9 FWHM is also cor- 
related with 7?NLR {rs = -0.494, Pr = 1.2 x 10"^). This 
leads to a a remarkably strong correlation of tnlr with 
Mbh {rs = 0.744, Pr = 3.8 x 10"*). This strong correlation 
stands in contrast with the lack of correlation of i?NLR with 
L/LEdd {rs = -0.002, Pr = 0.99). We note, however, that 
-Rnlr is not significantly correlated with Mbh {rs = 0.338, 
Pr = 3.3xl0~^). We do not have a simple plausible explana- 
tion for the strong J?nlr versus Mbh correlation. However, 
as mentioned above, one should note that this correlation 
is partly based on the /Jnlr versus L relation, which needs 
to be verified with higher quality spectra allowing complete 
detections. 



between U and rie, which produces a smaller than expected dis- 
persion in Uue, compared to the dispersion in each parameter. 
® The small spread in the BLR n~f can be understood in terms 
of the critical n-y required for dust sublimation (e.g. LD). 



4.5 The CF of the NLR 

Figure 8 presents a comparison of the single zone CF and 
the two zone CFin and CFout, for both the hard MF SED, 
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cates that most of the [O III] A5007 emission should be resolved 
with HST imaging. 
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Figure 7. The radial gradient in NLR density. The objects with 
upper limits may have a larger RnlRi and lower density, which 
would make this correlation stronger. 

and for the softer —1.5 slope SED (Section 4.1.3). In the 
two zone model with the soft SED we require the inner zone 
to produce a similar [Oiii] A4363/[Oin] A5007 ratio as in 
the hard SED model, which lead to log (7 = —1.4 for the 
inner zone (rather than log U = —1 with the hard SED). As 
already noted above (Section 4.1.3), the softer SED implies 
CF values which are higher by 1.74 compared to the hard 
SED. A similar effect is seen in Fig. 8 for the two zone so- 
lution, where CFin and CFout in the soft SED solution arc 
systematically higher than in the hard SED solution. For 
both SEDs CF and CFout are rather well correlated, but 
CFin can differ by a factor of 3-5 from CF. This is expected 



since CFin is largely set by [O iii] A4363, while both CF and 
CFout are largely set by [O iii] A5007. 

Independent constraints on the radial distribution of the 
covering factor of dusty gas in the NLR can be obtained from 
models for the IR SED of AGN. In the two zone model, the 
inner zone is within a factor of few of the dust sublimation 
radius (e.g. Fig. 8 in LD), and the associated dust radiation 
should thus peak at ~ 3 — 5 fxm, while the outer zone is at a 
^ 750 larger radius, and should thus peak at ^ 50—100 nm. 
The ratio CFout /CFin should then determine the shape of 
the far to near IR SED. Interestingly, Fig. 8 indicates that 
this ratio is systematically different for the soft and the hard 
SEDs, which suggests one may be able to obtain constraints 
on the ionizing SED based on the observed NLR emission 
and IR SED. The two zone model presented here is obviously 
highly simplified, and a more realistic model should invoke 
a continuous radial distribution of dusty gas. 

As (hscussed earlier, the large fraction of bolometric lu- 
minosity emitted in the IR (~ 0.3 — 0.5 in typical AGN), re- 
quires a correspondingly large CF of dusty gas (e.g. Sanders 
1989). Earlier studies noted that the CF deduced based on 
the EW of NLR non-resonant recombination lines is only a 
few percent, and concluded that there must exist an addi- 
tional NLR component with a large CF and U 3> 10^^ (Voit 
1992; LD, Netzer & Laor 1993). Our single zone model yields 
log CF~ -1 ± 0.5 (for the soft SED), which is still too low. 
A somewhat larger total CF is obtained with the two zone 
model, mainly due to our assumption of log (7 = — 1 in the 
inner zone. As briefly mentioned in Section 2, some of our 
solutions with log U < —2.2 will be degenerate with models 
with \ogU > — 1 (Fig. 3). Such high U models will require 
a higher CF, which may be in better agreement with the IR 
constraints on the CF. However, the best test for such high 
U models is through the various expected high ionization 
lines (originating in the surface layer where dust absorption 
is negligible). 

Finally, we note that the [O ill] A5007 line is one of the 
main NLR coolants over a rather wide range of Ue and U 
(e.g. Ferguson et al. 1997), and thus should be a good tracer 
for most of the NLR gas. The [O ill] A5007 Une becomes a 
negligible coolant at log U < —3.5 (e.g. Fig. 3), once most O 
is in O I and O II forms. However, a large covering factor of 
low U NLR gas without [O iii] A5007 emission can be ruled 
out in most objects based on the following argument. Such a 
low U component would account for a significant fraction of 
the H/3 line EW, and thus the inferred [O ill] A5007/H/3 ratio 
in the [O ill] A5007 emitting region will be significantly larger 
than the observed ratio. The commonly observed ratio of 
[Olll] A5007/H/3~ 10 is already quite close to the maximum 
possible value in photoionization models (e.g. Fig. 2, upper 
panel), and it does not allow more than ^ 1/2, of the H/3 
line to originate in non [O ill] A5007 emitting gas. 



5 CONCLUSIONS 

We use the [O iii] A5007 fine profile to measure the EW of 
the [O III] A4363 and H/3 emission lines in the BG92 sample 
of bright z < 0.5 AGN. All three lines are detected in 40 out 
of the 78 sample objects where optical spectra are available. 
The EW of the three lines are used, together with photoion- 
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Figure 8. Tlie covering factor of the outer (empty diamond) 
and inner (empty star) regions in the two-zones model versus the 
covering factor of the single-zone model. Upper panel, the relation 
for the assumed hard ionizing MF SED. Lower panel, the relation 
for the soft ionizing SED. A 1:1 relation is denoted by the dashed 
line in both panels. Note that that the ratio CFin/CFin in a 
given object, which sets the observed IR SED, is a function of the 
ionizing SED. 



ization models, to infer rie, U, and CF of the [Oiii] A5007 
emitting gas in the NLR. We find the main following results: 

1. The inferred and U in the single zone ap- 
proximation show a relatively small range (rie ~ 10^ — 
10®-^ cm-^ U ~ IQ-^-^ - IQ-*-^), which corresponds to the 
range whore the [O ill] A5007 cmissivity is maximized. 

2. The strength of the [O ill] A5007 line is mostly modu- 
lated by the range of CF of the photoionized gas. The range 
of [Oiii] A5007 emissivity per unit solid angle, set by rie, 
and [/, accounts for about half the range produced by the 
CF. Similar results are obtained for soft and hard ionizing 
SEDs, and for an extreme two zone approximation for the 
NLR. 

3. The NLR CF is inversely correlated with luminos- 
ity, consistent with some earlier suggestions, but it is not 
correlated with L/L^dd, unlike earlier suggestions. 

4. We find a rather tight radius luminosity relation, 
^?NLR = 4:0L%^^ pc, in the single zone approximation. 
We also find an unexpected strong correlation of i?NLR 



and Mbh- However, both relations may be induced by our 
[O III] A4363 detection limits . 

5. The observed line ratios can also be fit with a two 
zone model, where most of [O iii] A5007 originates in an 
He = 10'^ cm^'' component at i^^^LR = 7501/44*'* pc, and most 
of [O III] A4363 in an rie = 10^ cm~^ component at i?NLR = 
L44 pc. Significant constraints on the spatial distribution 
of the [O III] A5007 emission can be obtained through HST 
imaging of luminous Type I AGN, which should be able 
to resolve most of the [O III] A5007 emission if it occurs in 
Ue ~ 10^ cm~^ gas. 

Significantly higher quality optical spectroscopy of 
AGN, allowing complete detections of the [O ill] A4363 line, 
is required to establish the strength of the -Rnlr correlations 
with L and Mbii. 

Independent constraints on the radial distribution of the 
NLR gas covering factor can be deduced from the IR SED. 
These constraints, together with measurements of lower and 
higher ionization narrow lines, and more realistic photoion- 
ization models involving a radial distribution of dusty gas 
(e.g. Ferguson et al. 1997), can provide a better understand- 
ing of the origin of the observed extreme range in narrow line 
strengths in AGN. 
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also find that most of the [O ill] A5007 emission is signifi- 
cantly more compact, with half of the flux originating within 
-Rnlr = UOLlf^ pc, and typically a factor of two less in 
Type I AGN. This size is within a factor of 2-3 of our single 
zone J?NLR. 



NOTE ADDED IN PROOFS 

Schmitt et al. (2003) surveyed with HST the extended [Olll] 
emission in a sample of 60 nearby Type I/II Seyfert galax- 
ies. They find that the maximum detectable extent of the 
[O III] A5007 emission is given by -Rnlr = 790L44'*'* pc (as- 
summg a mean [O III] A5007 EW of 35 A to convert their 
[O III] A5007 luminosity to our continuum luminosity) , which 
is (fortuitously) very close to our expression for iiN^n,. They 



